Allosteric modulation of fluorescence revealed by hydrogen bond dynamics in a genetically encoded maltose biosensor.

Genetically encoded fluorescent biosensors (GEFBs) proved to be reliable tracers for many metabolites and cellular processes. In the simplest case, a fluorescent protein (FP) is genetically fused to a sensing protein which undergoes a conformational change upon ligand binding. This drives a rearrangement in the chromophore environment and changes the spectral properties of the FP. Structural determinants of successful biosensors are revealed only in hindsight when the crystal structures of both ligand-bound and ligand-free forms are available. This makes the development of new biosensors for desired analytes a long trial-and-error process. In the current study, we conducted µs-long all atom molecular dynamics (MD) simulations of a maltose biosensor in both the apo (dark) and holo (bright) forms. We performed detailed hydrogen bond occupancy analyses to shed light on the mechanism of ligand induced conformational change in the sensor protein and its allosteric effect on the chromophore environment. We find that two strong indicators for distinguishing bright and dark states of biosensors are due to substantial changes in hydrogen bond dynamics in the system and solvent accessibility of the chromophore.

Kinetic Barrier to Enzyme Inhibition Is Manipulated by Dynamical Local Interactions in E. coli DHFR.

Dihydrofolate reductase (DHFR) is an important drug target and a highly studied model protein for understanding enzyme dynamics. DHFR's crucial role in folate synthesis renders it an ideal candidate to understand protein function and protein evolution mechanisms. In this study, to understand how a newly proposed DHFR inhibitor, 4'-deoxy methyl trimethoprim (4'-DTMP), alters evolutionary trajectories, we studied interactions that lead to its superior performance over that of trimethoprim (TMP). To elucidate the inhibition mechanism of 4'-DTMP, we first confirmed, both computationally and experimentally, that the relative binding free energy cost for the mutation of TMP and 4'-DTMP is the same, pointing the origin of the characteristic differences to be kinetic rather than thermodynamic. We then employed an interaction-based analysis by focusing first on the active site and then on the whole enzyme. We confirmed that the polar modification in 4'-DTMP induces additional local interactions with the enzyme, particularly, the M20 loop. These changes are propagated to the whole enzyme as shifts in the hydrogen bond networks. To shed light on the allosteric interactions, we support our analysis with network-based community analysis and show that segmentation of the loop domain of inhibitor-bound DHFR must be avoided by a successful inhibitor.

Nonlocal Effects of Antibiotic-Resistance-Causing Mutations Reveal an Alternative Region for Targeting on FtsW-Penicillin-Binding Protein 3 Complex of Haemophilus influenzae.

Currently prescribed antibiotics target the catalytic sites of wild-type bacterial proteins; however, bacteria adopt mutations at this site, eventually leading to the emergence of resistance. Therefore, the identification of alternative drug binding sites is crucial, which requires knowledge of the dynamics of the mutant protein. Here, we set out to investigate the impact of a high-resistance-causing triple mutation (S385T + L389F + N526K) on the dynamics of a prioritized resistant pathogen, Haemophilus influenzae, using computational techniques. We studied penicillin-binding protein 3 (PBP3) and its complex with FtsW, which display resistance toward ß-lactam antibiotics. We showed that mutations displayed local and nonlocal effects. In terms of the former, the orientation of the ß-sheet, which surrounds the active site of PBP3, was impacted and the catalytic site was exposed to the periplasmic region. In addition, the flexibility of the ß3-ß4 loop, which modulates the catalysis of the enzyme, increased in the mutant FtsW-PBP3 complex. As for nonlocal effects, the dynamics of the pedestal domain (N-terminal periplasmic modulus (N-t)), i.e., the opening of the fork, was different between the wild-type and mutant enzymes. We showed the closed fork caused a greater number of residues to participate in the hypothesized allosteric communication network connecting N-t to the transpeptidase domain in the mutant enzyme. Finally, we demonstrated that the closed fork results in more favorable binding with ß-lactam antibiotics, particularly cefixime, suggesting that small therapeutics that can stabilize the closed fork of mutant PBP3 may lead to the development of more effective molecules to combat resistant bacteria.

Conformational multiplicity of bacterial ferric binding protein revealed by small angle x-ray scattering and molecular dynamics calculations.

This study combines molecular dynamics (MD) simulations with small angle x-ray scattering (SAXS) measurements to investigate the range of conformations that can be adopted by a pH/ionic strength (IS) sensitive protein and to quantify its distinct populations in solution. To explore how the conformational distribution of proteins may be modified in the environmental niches of biological media, we focus on the periplasmic ferric binding protein A (FbpA) from Haemophilus influenzae involved in the mechanism by which bacteria capture iron from higher organisms. We examine iron-binding/release mechanisms of FbpA in varying conditions simulating its biological environment. While we show that these changes fall within the detectable range for SAXS as evidenced by differences observed in the theoretical scattering patterns calculated from the crystal structure models of apo and holo forms, detection of conformational changes due to the point mutation D52A and changes in ionic strength (IS) from SAXS scattering profiles have been challenging. Here, to reach conclusions, statistical analyses with SAXS profiles and results from different techniques were combined in a complementary fashion. The SAXS data complemented by size exclusion chromatography point to multiple and/or alternative conformations at physiological IS, whereas they are well-explained by single crystallographic structures in low IS buffers. By fitting the SAXS data with unique conformations sampled by a series of MD simulations under conditions mimicking the buffers, we quantify the populations of the occupied substates. We also find that the D52A mutant that we predicted by coarse-grained computational modeling to allosterically control the iron binding site in FbpA, responds to the environmental changes in our experiments with conformational selection scenarios that differ from those of the wild type.

Inhibition of mutant RAS-RAF interaction by mimicking structural and dynamic properties of phosphorylated RAS.

Undruggability of RAS proteins has necessitated alternative strategies for the development of effective inhibitors. In this respect, phosphorylation has recently come into prominence as this reversible post-translational modification attenuates sensitivity of RAS towards RAF. As such, in this study, we set out to unveil the impact of phosphorylation on dynamics of HRASWT and aim to invoke similar behavior in HRASG12D mutant by means of small therapeutic molecules. To this end, we performed molecular dynamics (MD) simulations using phosphorylated HRAS and showed that phosphorylation of Y32 distorted Switch I, hence the RAS/RAF interface. Consequently, we targeted Switch I in HRASG12D by means of approved therapeutic molecules and showed that the ligands enabled detachment of Switch I from the nucleotide-binding pocket. Moreover, we demonstrated that displacement of Switch I from the nucleotide-binding pocket was energetically more favorable in the presence of the ligand. Importantly, we verified computational findings in vitro where HRASG12D/RAF interaction was prevented by the ligand in HEK293T cells that expressed HRASG12D mutant protein. Therefore, these findings suggest that targeting Switch I, hence making Y32 accessible might open up new avenues in future drug discovery strategies that target mutant RAS proteins.

Investigation of iron release from the N- and C-lobes of human serum transferrin by quantum chemical calculations.

Human serum transferrin binds ferric ions with high affinity and delivers them into cells via receptor-mediated endocytosis upon a decrease in pH in the endosome. Protonation events and conformational changes are known to play an important role in iron-release though the release is not yet fully understood. Human serum transferrin consists of two similar lobes which release iron at different rates. In this study, we investigate the iron binding sites of N- and C-lobes using quantum mechanical tools, particularly, the quantum chemical cluster approach. This study supports the inevitable role of axial tyrosine for the release of iron in quantum chemical models and provides valuable information about the proton transfer pathways for the protonation of Tyr188 and Tyr517 in N- and C-lobes, respectively. The calculations show that the release process is similar in both lobes; however, the energetic differences of the release process in N- and C-lobes, demonstrated for the first time, indicated that the release of iron in the N-lobe is thermodynamically favorable, in contrast to the one in the C-lobe.

Sibling rivalry among the ZBTB transcription factor family: homodimers versus heterodimers.

The BTB domain is an oligomerization domain found in over 300 proteins encoded in the human genome. In the family of BTB domain and zinc finger-containing (ZBTB) transcription factors, 49 members share the same protein architecture. The N-terminal BTB domain is structurally conserved among the family members and serves as the dimerization site, whereas the C-terminal zinc finger motifs mediate DNA binding. The available BTB domain structures from this family reveal a natural inclination for homodimerization. In this study, we investigated the potential for heterodimer formation in the cellular environment. We selected five BTB homodimers and four heterodimer structures. We performed cell-based binding assays with fluorescent protein-BTB domain fusions to assess dimer formation. We tested the binding of several BTB pairs, and we were able to confirm the heterodimeric physical interaction between the BTB domains of PATZ1 and PATZ2, previously reported only in an interactome mapping experiment. We also found this pair to be co-expressed in several immune system cell types. Finally, we used the available structures of BTB domain dimers and newly constructed models in extended molecular dynamics simulations (500 ns) to understand the energetic determinants of homo- and heterodimer formation. We conclude that heterodimer formation, although frequently described as less preferred than homodimers, is a possible mechanism to increase the combinatorial specificity of this transcription factor family.

DHFR Mutants Modulate Their Synchronized Dynamics with the Substrate by Shifting Hydrogen Bond Occupancies.

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize Cα-Cß bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-µs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

Computational strategies for protein conformational ensemble detection.

Protein function is constrained by the three-dimensional structure but is delineated by its dynamics. This framework must satisfy specificity of function along with adaptability to changing environments and evolvability under external constraints. The accessibility of the available regions of the energy landscape for a set of conditions and shifts in the populations upon their modulation have effects propagating across scales, from biomolecular interactions, to organisms, to populations. Developing the ability to detect and juggle protein conformations supplemented by a physics-based understanding has implications for not only in vivo problems but also for resistance impeding drug discovery and bionano-sensor design.

A trimethoprim derivative impedes antibiotic resistance evolution.

The antibiotic trimethoprim (TMP) is used to treat a variety of Escherichia coli infections, but its efficacy is limited by the rapid emergence of TMP-resistant bacteria. Previous laboratory evolution experiments have identified resistance-conferring mutations in the gene encoding the TMP target, bacterial dihydrofolate reductase (DHFR), in particular mutation L28R. Here, we show that 4'-desmethyltrimethoprim (4'-DTMP) inhibits both DHFR and its L28R variant, and selects against the emergence of TMP-resistant bacteria that carry the L28R mutation in laboratory experiments. Furthermore, antibiotic-sensitive E. coli populations acquire antibiotic resistance at a substantially slower rate when grown in the presence of 4'-DTMP than in the presence of TMP. We find that 4'-DTMP impedes evolution of resistance by selecting against resistant genotypes with the L28R mutation and diverting genetic trajectories to other resistance-conferring DHFR mutations with catalytic deficiencies. Our results demonstrate how a detailed characterization of resistance-conferring mutations in a target enzyme can help identify potential drugs against antibiotic-resistant bacteria, which may ultimately increase long-term efficacy of antimicrobial therapies by modulating evolutionary trajectories that lead to resistance.

Dynamic Community Composition Unravels Allosteric Communication in PDZ3.

The third domain of PSD-95 (PDZ3) is a model for investigating allosteric communication in protein and ligand interactions. While motifs contributing to its binding specificity have been scrutinized, a conformational dynamical basis is yet to be established. Despite the miniscule structural changes due to point mutants, the observed significant binding affinity differences have previously been assessed with a focus on two α-helices located at the binding groove (α2) and the C-terminus (α3). Here, we employ a new computational approach to develop a generalized view on the molecular basis of PDZ3 binding selectivity and interaction communication for a set of point mutants of the protein (G330T, H372A, G330T-H372A) and its ligand (CRIPT, named L1, and its T-2F variant, L2) along with the wild type (WT). To analyze the dynamical aspects hidden in the conformations that are produced by molecular dynamics simulations, we utilize variations in community composition calculated based on the betweenness centrality measure from graph theory. We find that the highly charged N-terminus, which is located far from the ligand, has the propensity to share the same community with the ligand in the biologically functional complexes, indicating a distal segment might mediate the binding dynamics. N- and C-termini of PDZ3 share communities, and α3 acts as a hub for the whole protein by sustaining the communication with all structural segments, albeit being a trait not unique to the functional complexes. Moreover, α2 which lines the binding cavity frequently parts communities with the ligand and is not a controller of the binding but is rather a slave to the overall dynamics coordinated by the N-terminus. Thus, ligand binding fate in PDZ3 is traced to the population of community compositions extracted from dynamics despite the lack of significant conformational changes.

N-Terminus of the Third PDZ Domain of PSD-95 Orchestrates Allosteric Communication for Selective Ligand Binding.

PDZ domains constitute common models to study single-domain allostery without significant structural changes. The third PDZ domain of PSD-95 (PDZ3) is known to have selective structural features that confer unique modulatory roles to this unit. In this model system, two residues, H372 directly connected to the binding site and G330 holding an off-binding-site position, were designated to assess the effect of mutations on binding selectivity. It has been observed that the H372A and G330T-H372A mutations change ligand preferences from class I (T/S amino acid at position -2 of the ligand) to class II (hydrophobic amino acid at the same position). Alternatively, the G330T single mutation leads to the recognition of both ligand classes. We have performed a series of molecular dynamics (MD) simulations for wild-type, H372A, and G330T single mutants and a double mutant of PDZ3 in the absence and presence of both types of ligands. With the combination of free-energy difference calculations and a detailed analysis of MD trajectories, "class switching" and "class bridging" behavior of PDZ3 mutants, as well as their effects on ligand selection and binding affinities are explained. We show that the dynamics of the charged N-terminus plays a fundamental role in determining the binding preferences in PDZ3 by altering the electrostatic energy. These findings are corroborated by simulations on N-terminus-truncated versions of these systems. The dynamical allostery orchestrated by the N-terminus offers a fresh perspective to the study of communication pathways in proteins.

A Coarse-Grained Methodology Identifies Intrinsic Mechanisms That Dissociate Interacting Protein Pairs.

We address the problem of triggering dissociation events between proteins that have formed a complex. We have collected a set of 25 non-redundant, functionally diverse protein complexes having high-resolution three-dimensional structures in both the unbound and bound forms. We unify elastic network models with perturbation response scanning (PRS) methodology as an efficient approach for predicting residues that have the propensity to trigger dissociation of an interacting protein pair, using the three-dimensional structures of the bound and unbound proteins as input. PRS reveals that while for a group of protein pairs, residues involved in the conformational shifts are confined to regions with large motions, there are others where they originate from parts of the protein unaffected structurally by binding. Strikingly, only a few of the complexes have interface residues responsible for dissociation. We find two main modes of response: In one mode, remote control of dissociation in which disruption of the electrostatic potential distribution along protein surfaces play the major role; in the alternative mode, mechanical control of dissociation by remote residues prevail. In the former, dissociation is triggered by changes in the local environment of the protein, e.g., pH or ionic strength, while in the latter, specific perturbations arriving at the controlling residues, e.g., via binding to a third interacting partner is required for decomplexation. We resolve the observations by relying on an electromechanical coupling model which reduces to the usual elastic network result in the limit of the lack of coupling. We validate the approach by illustrating the biological significance of top residues selected by PRS on select cases where we show that the residues whose perturbation leads to the observed conformational changes correspond to either functionally important or highly conserved residues in the complex.

Structural analysis of the PATZ1 BTB domain homodimer.

PATZ1 is a ubiquitously expressed transcriptional repressor belonging to the ZBTB family that is functionally expressed in T lymphocytes. PATZ1 targets the CD8 gene in lymphocyte development and interacts with the p53 protein to control genes that are important in proliferation and in the DNA-damage response. PATZ1 exerts its activity through an N-terminal BTB domain that mediates dimerization and co-repressor interactions and a C-terminal zinc-finger motif-containing domain that mediates DNA binding. Here, the crystal structures of the murine and zebrafish PATZ1 BTB domains are reported at 2.3 and 1.8 Šresolution, respectively. The structures revealed that the PATZ1 BTB domain forms a stable homodimer with a lateral surface groove, as in other ZBTB structures. Analysis of the lateral groove revealed a large acidic patch in this region, which contrasts with the previously resolved basic co-repressor binding interface of BCL6. A large 30-amino-acid glycine- and alanine-rich central loop, which is unique to mammalian PATZ1 amongst all ZBTB proteins, could not be resolved, probably owing to its flexibility. Molecular-dynamics simulations suggest a contribution of this loop to modulation of the mammalian BTB dimerization interface.

Perturb-Scan-Pull: A Novel Method Facilitating Conformational Transitions in Proteins.

Conformational transitions in proteins facilitate precise physiological functions. Therefore, it is crucial to understand the mechanisms underlying these processes to modulate protein function. Yet, studying structural and dynamical properties of proteins is notoriously challenging due to the complexity of the underlying potential energy surfaces (PES). We have previously developed the perturbation-response scanning (PRS) method to identify key residues that participate in the communication network responsible for specific conformational transitions. PRS is based on a residue-by-residue scan of the protein to determine the subset of residues/forces which provide the closest conformational change leading to a target conformational state, inasmuch as linear response theory applies to these motions. Here, we develop a novel method to further evaluate if conformational transitions may be triggered on the PES. We aim to study functionally relevant conformational transitions in proteins by using results obtained from PRS and feeding them as inputs to steered molecular dynamics simulations. The success and the transferability of the method are evaluated on three protein systems having different complexities of motion on the PES: calmodulin, adenylate kinase, and bacterial ferric binding protein. We find that the method captures the target conformation, while providing key residues and the optimum paths with relatively low free energy profiles.

High-Order Epistasis in Catalytic Power of Dihydrofolate Reductase Gives Rise to a Rugged Fitness Landscape in the Presence of Trimethoprim Selection.

Evolutionary fitness landscapes of several antibiotic target proteins have been comprehensively mapped showing strong high-order epistasis between mutations, but understanding these effects at the biochemical and structural levels remained open. Here, we carried out an extensive experimental and computational study to quantitatively understand the evolutionary dynamics of Escherichia coli dihydrofolate reductase (DHFR) enzyme in the presence of trimethoprim-induced selection. To facilitate this, we developed a new in vitro assay for rapidly characterizing DHFR steady-state kinetics. Biochemical and structural characterization of resistance-conferring mutations targeting a total of ten residues spanning the substrate binding pocket of DHFR revealed distinct changes in the catalytic efficiencies of mutated DHFR enzymes. Next, we measured biochemical parameters (Km, Ki, and kcat) for a mutant library carrying all possible combinations of six resistance-conferring DHFR mutations and quantified epistatic interactions between them. We found that the high-order epistasis in catalytic power of DHFR (kcat and Km) creates a rugged fitness landscape under trimethoprim selection. Taken together, our data provide a concrete illustration of how epistatic coupling at the level of biochemical parameters can give rise to complex fitness landscapes, and suggest new strategies for developing mutant specific inhibitors.

Tuning Interaction Parameters of Thermoplastic Polyurethanes in a Binary Solvent To Achieve Precise Control over Microphase Separation.

Thermoplastic polyurethanes (TPUs) are designed using a large variety of basic building blocks but are only synthesized in a limited number of solvent systems. Understanding the behavior of the copolymers in a selected solvent system is of particular interest to tune the intricate balance of microphase separation/mixing, which is the key mechanism behind the structure formation in TPUs. Here, we present a computationally efficient approach for selecting TPU building blocks and solvents based on their Flory-Huggins interaction parameters for a precise control over the microphase separation/mixing. We first cluster eight soft segments (PEO, PPO, PTMO, PBA, PCL, PDMS, PIB, or PEB) used frequently in TPUs into three categories according to the strength of their interactions with the binary solvent THF/DMF. We then perform a comprehensive set of dissipative particle dynamics simulations of the TPUs in a range of solvent ratios. This enables us to demonstrate the emergence of the unusual channel-like structures in a narrow range of parameters and to determine the critical interactions operative for obtaining either microphase separated or mixed structures. The findings are supported by thermodynamic arguments. The approach developed here is useful for designing novel TPUs with well-defined conformational characteristics, controlled morphologies, and advanced functional properties.

Employment of Iron-Binding Protein from Haemophilus influenzae in Functional Nanopipettes for Iron Monitoring.

Because of the serious neurologic consequences of iron deficiency and iron excess in the brain, interest in the iron status of the central nervous system has increased significantly in the past decade. While iron plays an important role in many physiological processes, its accumulation may lead to diseases such as Huntington's, Parkinson's, and Alzheimer's. Therefore, it is important to develop methodologies that can monitor the presence of iron in a selective and sensitive manner. In this paper, we first showed the synthesis and characterization of the iron-binding protein (FBP) from Haemophilus influenzae, specific for ferrous ions. Subsequently, we employed this protein in our nanopipette platform and utilized it in functionalized nanoprobes to monitor the presence of ferrous ions. A suite of characterization techniques: absorbance spectroscopy, dynamic light scattering, and small-angle X-ray scattering were used for FBP. The functionalized Fe-nanoprobe calibrated in ferrous chloride enabled detection from 0.05 to 10 µM, and the specificity of the modified iron probe was evaluated by using various metal ion solutions.

Computational Methods for Efficient Sampling of Protein Landscapes and Disclosing Allosteric Regions.

Methods developed toward computational exploration of protein landscapes have become standardized tools to assess biophysical experimental findings. They are also used on their own right to discover the workings of the protein as a molecular machine, potential sites of interest for protein functioning, allosteric regions in proteins, and communication pathways between different sites on a protein. With the development of reliable force fields that describe interactions in biomolecules, molecular dynamics (MD) simulations have become the prime tool for this purpose. While it is now straightforward to carry out MD simulations up to microseconds with current computers readily available to researchers, many processes of biological interest occur on several of orders of magnitudes slower timescales. Thus, the latter problems are attackable through MD by a handful of researchers that have access to the most powerful computers. Alternatively, physics-based methods to interrogate the protein energy landscape are in continuous development to circumvent this problem. In addition to opening the routes for advancement to a large number of researchers that have access to modest computational resources, they have the advantage of providing an understanding of the mechanisms that govern protein dynamics. Here we discuss network-based approaches geared toward understanding protein dynamics. These include (i) construction of residue networks which view proteins as networks of nodes connected through local interactions and (ii) construction of proteins as elastic networks whose modes of motion may be manipulated to achieve allowed conformational changes. Limitations of the methods as well as opportunities for future exploitation are described.

MODE-TASK: large-scale protein motion tools.

Summary: MODE-TASK, a novel and versatile software suite, comprises Principal Component Analysis, Multidimensional Scaling, and t-Distributed Stochastic Neighbor Embedding techniques using Molecular Dynamics trajectories. MODE-TASK also includes a Normal Mode Analysis tool based on Anisotropic Network Model so as to provide a variety of ways to analyse and compare large-scale motions of protein complexes for which long MD simulations are prohibitive. Beside the command line function, a GUI has been developed as a PyMOL plugin. Availability and implementation: MODE-TASK is open source, and available for download from https://github.com/RUBi-ZA/MODE-TASK. It is implemented in Python and C++. It is compatible with Python 2.x and Python 3.x and can be installed by Conda. Supplementary information: Supplementary data are available at Bioinformatics online.